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Abstract 



John Holte |16| introduced a family of "amazing matrices" which give the transition 
probabihties of "carries" when adding a hst of numbers. It was subsequently shown 
that these same matrices arise in the combinatorics of the Veronese embedding of com- 
mutative algebra [H [5J [7] and in the analysis of riffle shuffling [51 [7] . We find that the 
left eigenvectors of these matrices form the Foulkes character table of the symmetric 
group and the right eigenvectors are the Eulerian idempotents introduced by Loday 
|20j in work on Hochschild homology. The connections give new closed formulae for 
Foulkes characters and allow explicit computation of natural correlation functions in 
the original carries problem. 

Keywords: Foulkes character, carry, Eulerian idempotent, symmetric group 
AMS 2010 subject classifications: primary 20C30; secondary 60C05, 60J10. 

1 Introduction 

When n numbers are added in the usual way, "carries" accrue along the way. For example, 
working base h = 10, the display shows the carries along the top when n = 3 ten-digit 
numbers are added: 



Here the carries (reading right to left in bold print) are kq = 0, ki = 2, K2 = 1, K3 = 2, . . . . 
When n numbers are added, the carries can be 0, 1, 2, . . . , n — 1. If the digits are chosen 
uniformly at random in {0, 1, . . . , 6 — 1}, it is easy to see that the carries form a Markov 
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chain: the chance that the next carry is j given the past carries only depends on the last 
carry. Thus the distribution of carries is determined by the transition matrix 



M{i,j) = chance {next carry is j \ last carry is i}. 
The carries process was studied by Holte [16] who showed 

/ , , , 1^"^^^, ,,/n + l\ /n-l-i + (j + l-l)b\ 

(1.1) M{^,J) = - {-iy( ^ M, 0<^,J<n-l. 

For example, when n = 3, the matrix is 



462 


-4 


62 - 36 + 2 


462 


+ 2 


62-1 


462 


-4 


62 + 36 + 2 



1 

662 



Holte found the eigenvalues, eigenvectors, and many amazing properties of these matrices. 

Work of [H ini [7] shows that the same matrix arises in the analysis of the Gilbert- 
Shannon-Reeds method of shuffling cards and in the Hilbert series of the Veronese embed- 
ding of projective varieties. 

The main results of the present paper identify a different area where the matrix appears. 
The left eigenvectors of the matrix are the Foulkes characters of the symmetric group. The 
right eigenvectors are the Eulerian idempotents that occur in the study of free Lie algebras 
and Hochschild homology. We obtain new closed-form expressions for these characters. 



Section 2 gives background on Foulkes characters and presents some new results for left 



eigenvectors. Section 3 does the same for the right eigenvectors and applies some of the new 
formulae to the original carries process, giving the variance and covariance of the number of 
carries. [Section 4 gives another connection between representation theory of the symmetric 
group (the RSK correspondence) and carries. 



2 Foulkes characters 



This section introduces the Foulkes characters of the symmetric group and some of their 
properties (Section 2.1 ). It shows that the Foulkes characters are the left eigenvectors of the 



transition matrix M of (1.1) (Section 2.2). This connection is used to prove a branching 
rule (from 5„ to 5„_i) and a closed- form formula for Foulkes characters (Section 2.3). 



2.1 Background on Foulkes characters 



Foulkes characters were discovered by Foulkes ^0\ as part of the study of the descent 



patterns in the permutation group. They are developed in [18] and [T7] gives a readable 



textbook treatment. Gessel and Reutenauer |14| use them to enumerate permutations by 
descents and conjugacy classes; see [5] for a probabilistic interpretation of these results. 



Stanley [26] uses Foulkes characters to develop enumerative results for alternating permu- 



tations by cycle type. 

Recall that a permutation a £ Sn has a descent at i if a{i + 1) < cr(i). The set of places 
where descents occur is D(a) C [n — 1]. For example, if a = 4512376, D{a) = {2,6}. If 
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U C [n — 1] is any set, Foulkes suggested constructing a ribbon shape (also called a rim 
hook) R{U) beginning with a single box and sequentially adding the next box below the 
last box if i G [7, and to the left of the last box if i ^ U, 1 < i < n — 1. Thus, if f7 = {2, 6}, 
boxes are built up as follows: 



™ F cff off cuff gzff 

The final skew shape will have n boxes and be the lower rim of a partition a; in this example, 
a is 5,4,1, and the ribbon shape is 5,4, 1\3: 




Labeling the boxes in the ribbon shape by all ways they can be sequentially removed from 
a and reading this from right to left and top to bottom gives all permutations with the 
original U as descent set. For example, removing boxes in the order shown as 









4 




7 


6 


5 


3 
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gives 1435672. The skew shape R{U) corresponding to ?7 C [n — 1] gives a skew character 
Ii(JJ^ = a\l3 and is an irreducible character of the symmetric group Sn, the 
coefficient of ™ x^^^^ is {x^ ' X^\x°') (see |21l Sect. 1.7]). From the development above, 
the dimension of x'^^^'^ is the number of permutations with descent set U. Solomon [251 
Sect. 6] describes a related construction of MacMahon in his work on Simon Newcomb's 
problem. 

For fixed k, < k < n — 1, the Foulkes character x"'*" is defined as the sum of x^^^'^ 
over all U with n — k — 1 descents. It follows that the dimension of x^'^ is the Eulerian 
number A[n, k), the number of permutations with k descents. Foulkes showed that x^'^^i'^) 
only depends on a through the number of cycles in a. In particular, 

(2.1) x"'"""*^ is the trivial character; 

(2.2) x^'^ is the alternating character. 

Most importantly, letting xj''^ denote the value of the Foulkes character on permutations 
with j cycles (so the dimension Xn'^ = A{n, k)), 

/n n\ "i^ n— l,fc— 1 n—l,k -r ■ ^ 7 ^ n 

(2.3) Xj =Xj -Xj ifj<n, fc>0. 

This, with the starting value Xi '^ = li gives an efficient way to build a Foulkes character 
table. Let k = 0, 1, . . . , n — 1 index the rows and j = n,n — 1, . . . ,1 index the columns. 



Table 1 gives the example when n = 5. 
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2 


1 





1 


-1 


1 


-1 


1 


1 


26 


-10 


2 


2 


-4 


2 


66 





-6 





6 


3 


26 


10 


2 


-2 


-4 


4 


1 


1 


1 


1 


1 



Table 1: Foulkes character table for n = 5. 



Further properties of Foulkes characters appear in Kerber and Thiirlings |18j : 

(2.4) > ^ Ai < A: + 1, a; < n - A:. 

(2.5) (x«'^x^■+'•'""")>0«^j = A;. 

(2.6) The x"''^ are Unearly independent. 

(2.7) If X : 'S'n — IS denotes a character, depending only on the number of cycles, then 



X 



E 



X 



(id) 



X ' ■ 



Thus the x" form a Q basis for the characters that only depends on the number of cycles. 

Hidden in the proof of (2.7): the hook character x*"*"^'^" ' is the only hook oc curing in 

rn-l\ ^ „„i„i„j f„„i „ ;„ ISqIq^q^ 



X"'* and it occurs with multiplicity its degree (J'' ■ ). A related fact appears in 



Th. 4]. [Kerber and Thiirlings 



acting on [M]": 
(2.8) 

has the decomposition 



further determine the permutation character for Sn 
Xm{(7) ■■=M* '=^=1'^^'° ('^^ 



XM 



E 

k 



n 



The x"''^ cire usually not irreducible, and [Kerber and Thiirlings |18| give an interesting 
combinatorial rule for decomposing x^*"^^ (and thus x"''')- They show that x"'*^ is a sum 
of Lefschetz characters for Sn acting on 21"]. 

Marty Isaacs conjectured that n\ divides the determinant of the Foulkes character table. 
In fact the following is true: 



(2.9) The determinant of the n x n matrix with k,j entry x"''^ is n!(n 



l)!---2!. 



Proof. Construct an n x n matrix A from the (n — 1) x (n — 1) Foulkes character table by 
adding a left column consisting of the partial sums of the Eulerian numbers A{n, 0),A{n, 0) + 
A(n, 1), . . . , n! and filling out the rest of the top row with zeros. Thus, when n = 5, 



A 



/120 











o\ 


119 


1 


-1 


1 


-1 


93 


11 


-3 


-1 


3 


27 


11 


3 


-1 


-3 


V 1 


1 


1 


1 


1 / 
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The 4x4 matrix in the lower right corner is the Foulkes character table for n = 4. The first 
column entries are the partial sums of the Eulerian numbers 1, 26, 66, 26, 1. In particular, 
the (1,1) entry is n!, so by induction the determinant of A is n\{n — 1)! • • • 2!. 

The n X n Foulkes character table is constructed from A as follows: in A, subtract row 



2 from row 1, then row 3 from row 2, and so on. The recurrence (2.3) shows this gives the 
n X n Foulkes character table. □ 

Marty Isaacs observes that while the {x"''^} i^ot disjoint, they sum up to the regular 
character of S'„; see [231 Th. 2] for a proof. Alas, this does not seem to be enough to have 
the nice theory of supercharacters [1] carry over, but the parallels are intriguing. Further 



properties of Foulkes characters are given in Section 2.3 after the connection with the carries 
transition matrix is developed. 

Rim hook characters are a basic construction of representation theory of Sn', see [3l I19| 
and the references there. They are also available for other Coxeter groups |23j . Foulkes' 
innovation, showing that sums of these characters have interesting properties, has not been 
explored for general type. 

2.2 The connection with carries 

We noticed from Holte's paper |16) that the matrix of left (row) eigenvectors for the carries 
Markov chain on {0, 1, 2, 3, 4} (e.g., working base 5) is 



n 


26 


66 


26 


1\ 


1 


10 





-10 


-1 


1 


2 


-6 


2 


1 


1 


-2 





2 


-1 


V 


-4 


6 


-4 


1/ 



Comparing this with Table 1 , the Foulkes character table, leads to the following result. 



Theorem 2.1. Let vf'j denote the jth entry of the ith left eigenvector of the carries matrix 
for addition of n numbers base b (here < i,j < n — 1, and the eigenvalues are 1/b^)- Then 



n,n—j—l 
n—i 



Proof. The first case is that i = 0. From |T6], Vqj = A{n,j). From the dimension formula, 

n,n—j—l 



Xn " = A{n,n—l—j). By symmetry of the Eulerian numbers, ^(n,j) = ^(n, n — 1— j), 
so the theorem follows in the first case. 

n,0 



The second case is that j 



n 



1. By equation (2.2), x,' = (—1)" Thus we need 



to show that f , 



n—j,n—l 



(—1)" K By Holte's formula for the left eigenvectors of the carries 



chain [161 P- 143], it follows that 



V, 



n-j,n-l 



Y^i-lY(^+'']in-ry 
r=0 ^ 
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The result now follows by induction, since 



J2{-lY(''^/)in-ry-\n-r) 



n 

r=0 ^ ^ r=0 

_ ^ 1) " (n - r)^-i 

r=l ^ ^ 

n-2 / V 

+ (n + 1) J^(-l)" ( j (n - 1 - ry-^ 

^_l)n-i+l + (^+l)(_l)n-, 
(-1)"-^'. 



For the remaining cases, i > and j < n — 1. By the recursive formula (2.3), it is 
enough to show that 



,,n _ n-1 _ n-1 

— '^1-1,3 ^i-l,j-l 



for i > 0, j < n — 1. From page 144 of |16j . 

(I \ n—i 
x—\ {l-x)-\ 



Clearly 

(7 \ n—i 
x—j (l-x)-^ 

(J \ n—i 
x—\ {l-x)-^ 

\ n—i 

X—j {i-xy^ 

which implies the result. □ 
2.3 Some consequences 

In [T6j, Holte gave a closed formula for the left eigenfunctions, 
(2.10) ^«. = ^(-lfr+ + 

r=0 V / 

Thus we get an apparently new formula for the Foulkes characters. 
Corollary 2.2. 

n-k , ^ 

I W 

V 



xf = Y^{-lY('''^'){n-k-ry 



In rereading Foulkes |10^ Sect. 4] we found the formula (Th. 4.1), 



n,k 



n-j 



E(-ir 



r=0 



n-j 



A{j, k + j + r — n). 



This seems a httle less direct than Corollary 2.2 Corollary 2.2 gives a direct proof of the 
following restriction formula of Foulkes characters from Sn to Sn~i- 

Corollary 2.3 ([10, Cor. 4.6]). 



n—l.k 



n-l,k-l 



< < n - 1. 



Remarks We first learned Corollary |2.3| from Marty Isaacs, who both observed it and 
showed that it follows from Corollary |2.2[ The well-known recursion formula for the Eulerian 
numbers A{n, k) = (k + l)A{n — l,k) + {n — k)A(n — 1, A; — 1) is the special case of evaluation 
at the identity. Thus Corollary |2.3| represents a "categorification" of this recurrence. 



Proof. The required formula translates to 

Xf = {k + l)xp,''' + {n- k)xY-'i'-\ < < n - 1, i > 1. 



From Corollary 2.2 



(k + Mz'i' + {n- k)xW'~' 



3- 

n~k~l , s 

(A; + l) (-ir(^)(n-A;-r-l)^-i + (n-fc)xp^''' 

n—k , \ 

(fe + l)^(-l)'^-M/; (n-A;-r)^-i 

r=l ^ ^ 



n—k , \ 

+ (n - kY + (n - A;) J]](-l)'' y\{n-k- -V^^ 

n—k / \ 

(n-A;)J +^(-1)^-1 f ^ \{n-k-ry-^ 

r=l ^ ' 

n—k / \ 



r 



11 — r<, y \ 

r=l ^ ^ 

n—k / \ n—k / \ 

+ n ^(-1)^' \^\{n-k- ry-^ - kY^{-\Y y\{n - k - ry-^ 
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fL— K / \ 

, - ky + n ^(-1)'' r]{n-k-ry 

r=l ^ ^ 

n—k / s 



+ ^ j(n- fe-ry-i ^ ) (n - A; - r^-^ 

r=l V / V 



n-fe , ^ 

^(-irr+')(n-^-rri 

r=l ^ 



in — k — ry 



(n-fcv+5:(-irr; 

r=l ^ 

n— fc / \ n—k / \ 

■^nY^{-\Y\^\(n-k-ry-^ ^Y^^-Xy-^i " J (n - A: - r)^-^ 

r-=l ^''^ r=l ^ 



n—k 



Y^{-\y[n -r)y'^^]{n-k- ry-^ 

r=l ^ 



n 
r-1 



[n — r] 



n + 1 



n,k 

where the final equahty used the identity 



n 



, , n + 1 , 
n-r 1=0. 
r — 1 / \ r 



□ 



3 Riffle shuffle idempotents and right eigenvectors of the car- 
ries matrix 

This section simplifies Holte's formula for the right eigenvectors of the carries matrix, and 



relates these eigenvectors to representation theory of the symmetric group (Section 3.1). 



group ( 


Section 3.1 


Section 3.2 



The eigenfunctions are used to compute basic things about carries in Section 3.2 
3.1 Right Eigenfunctions 

To begin, let denote the value of the jth right eigenvector of the carries chain (0 < 

j < n — 1, eigenvalues l/b^) evaluated at i {0 < i < n — 1). We also let s{n, k) be the Stirling 
number of the first kind, defined as (—1)""'^ multiplied by the number of permutations on 
n symbols with k cycles. It can also be defined by the equation 



(3.1) 



x{x — 1) . . . (x — n + 1) = s(n, k)x^ . 



k>0 
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Theorem 4 of Holte [T6| shows that 



k=n—j 



k-in-j) 



where is taken to be 1. Note that (i) is a polynomial in i of degree j. For example, 
<(i) = l, <(i)=n(n-l-i)-(^). 

The next theorem gives a simpler formula for u^{i). 



3 

Theorem 3.1. 



fc>0 

n! • Coefficient of x""-' in ' ~'~ 

n 



Proof. By Theorem 4 of Holte |16j . is n\ multiplied by the the i, j entry of the inverse 

of the matrix of left row eigenvectors. From this and equation (2.10), proving the first 
equality of the theorem is equivalent to proving that 



Z>0 fc>0 ■ ^ ^ r=0 ^ ^ 



Now 



"'EE=^("„-:r)B-<r)«--) 



!E("„-:r)E(-r("r)^E*"-o(i + i-o 

'E("„-:r)B-K"r)^t»(M)o-..-.)' 



«'e("::"/)»-'''^"^'^^'^'"' 

A;>0 ^ ^ r>0 

+ 1\ v-^ /n — z — iX/i + l — r 
n! >(-!)■ I I > I 

n — K 



'E(-i)'- E 

r>0 ^ ' k>0 

!E(-ir("r)("^^"""'') 



r>0 

n\6ij. 



The third equality used equation (3.1 ), the fifth equality used the basic identity (") 
^a+b^ PI p. 12], and the final equality is from page 147 of P^. 
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To prove the second equality of the theorem, write 



X + n — i — l\ /x\/n — i — 1 

n J ^-^ \k J \ n — k 



and use equation (3.1). □ 



Remark Let k be elements of the symmetric group algebra defined by the equation 



kj^ V- (x + n-d{w) 

X £jn_h = > 



fc=l W^Sn 



n 



where d{'w) denotes the number of descents of w. By work of Garsia and Reutenauer 
|12j . these are orthogonal idempotents of the symmetric group algebra whose sum is the 
identity. They also arise in the theory of riffle shuffling |2j and in Hochschild homology [15j . 
Their images under the sign map are known as Eulerian idempotents. In this version, they 
were discovered by Gerstenhaber-Schack [13] to give Hodge decompositions of Hochschild 
homology. They have been developed by Loday [20] for cyclic homology. Patras |22j gives 
an unusual treatment involving decompositions of the n-cube into simplices. For a textbook 
treatment, see Weibel [271 Sect. 9.4.3]. 

The eigenvectors of the carries and descent matrix lift to eigenvectors of the full riffle 
shuffle matrix. These in turn are identified in Denham [^ and Diaconis-Ram p]. 

Clearly the value of -En.fc on a permutation depends only on its number of descents. 
Letting E^^kid) denote the value of En^k on a permutation with d descents, we have the 



following corollary of Theorem 3.1 
Corollary 3.2. 

Uj{i) = n\En,n-j{i)- 
It would be nice to have a more conceptual proof of Corollary |3.2 



3.2 Applications 

This section gives some applications of the explicit form of the right eigenvectors of the 
carries chain for the addition of n numbers base b. We note that another application (to 
lower bounding the convergence rate of the carries chain) appears in [7]. 

The transition matrix of the carries chain is viewed as a linear operator on functions in 
the usual way: K f[x\ = Yl,y K{x, y)f{y)- Let Kr denote the value of the carry from column 
r — 1 to column r. 

Proposition 3.3. Suppose that the carries chain is started from its stationary distribution 
IT. Then for n>2, the covariance 

^ . . 1 n + 1 

Cov(Ko«^r) = ^—[^ > 0- 

Proof. Clearly 

E^KQKr) = i7r(«)£'(Kr|Ko = i) ■ 
i 
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By Theorem 3.1 , u" is a right eigenvector of the carries chain with eigenvalue 1/6. Dividing 
by — 1/n, we have the eigenvector f{i) = i — 



E(Kr\KQ 



n 



1 



1 



n 



1\ n 

+ 



1 



It follows that 
(3.2) 



E{K,QKr) = 



n 



+ 



n 



7r(i). 



and 



n+l 
12 ■ 



Prom [TB] or [7], the stationary distribution of the carries chain is 7r{i) = A{n,i)/n\ 
(here A(n,i) is the number of permutations on n symbols with i descents), and it is well 
known that for n > 2 the mean and variance of the Eulerian numbers are 
This, together with equation (3.2), gives 

E{KoKr) = — — tt; h 



1 

¥ 



12 
1 n + 1 



n 



1 

If 



n 



+ 



n 



12 



+ 



n 



□ 



and the result follows since E{Ko)E{Kr) 

A similar calculation allows us to compute the covariance started from the state 0: 
Proposition 3.4. Suppose that the carries chain is started from 0, and that n>2. Then 

1. E{n,,)={l-^)^. 



= (i-^) 



1 ^ n+l 
W) 12 • 



2. Var(Kr) = (1 

3. CoV^KgjKg^rJ — W^[2^ 

Proof. Parts 1 and 2 are proved in Theorem 4.1 of [6]. For part 3, arguing as in the proof 



1 (1 - 5^) > 0. 



of Proposition 3.3 



E{KsKs+r) = ^iP(Ks = i)E{Ks+r\Ks = i) 



n 



1 \ n 

+ 



1 



2 

n — 1 



E[k,]. 



The result follows from parts 1 and 2 after elementary simplifications. 



□ 



The eigenvectors of the carries matrix can be used to give a simple proof that the sum of 
^n,k giygg ^]^g regular character. Holte [16] shows that if V is the matrix whose rows are left 
eigenvectors and U is the matrix whose columns are right eigenvectors, then VU = n\ x I. 
Just looking at products involving the first column of U (which is identically one) we get 
that the x"''^ sum up to the regular character. 
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It is instructive to see the problems encountered in trying to use the available eigen- 
structure to bound the rate of convergence of the carries chain to its stationary distribution 
A(n, j)/n\. Let Mb{i,j) be the transition matrix (1.1) corresponding to adding n numbers 
base b. Let M^{i,j) be the kth power of this matrix. From elementary linear algebra, 



AinJ) 



n—l ^ 



a=l 



where la-,ra are the left and right eigenvectors of Mi, normed so that /qO) 
ro(j) = 1) for < J < n — 1. Here 



(and 



1 



= —,Xn-a^ ^ with 



n,n—j—l 



the Foulkes characters of 



Section 2 



^a(j) = """(j) with u^{j) given in Theorem 3.1 



The total variation distance to stationarity, starting at i = 0, after k steps is 

n-l 

2 



^ n— 1 



3=0 



A{n,j] 



From the formulae above, 



A{n,j) 



nl 



n— 1 ^ 



a=l 



While bounding this is feasible, it is a bit of a mess. In [71 Sect. 3], a different representation 
is used to prove that the carries chain is close to stationarity after order | log^(n) steps. 



4 Carries and the RSK correspondence 

In this section we use the RSK correspondence to derive a generating function for descents 
after a 5''-riffle shuffle on a deck of n cards. See for background on riffle shuffles. By 
a main result of [6j, this is equal to the generating function for the rth carry Kr when n 



numbers are added base b (and one can give another proof of Theorem 4.1 using Holte's 
formula for P{Kr = i)). 

Theorem 4.1. Let w be produced by a b^' -riffle shuffle on a deck of n cards, and let d{w) 
denote the number of descents of w. Then 

,4.1) ('-'=--'). 

W&Sn k>l ^ ^ 

Proof. Let w he a permutation produced by a 5'" riffle shuffle. The RSK correspondence 
associates to w; a pair of standard Young tableaux {P{w),Q{w)) of the same shape. More- 
over, there is a notion of descent set for standard Young tableaux, and by Lemma 7.23.1 
of [25], the descent set of w is equal to the descent set of Q{w). It is known from [TT] 
that if w is produced by a h'' shuffle, then the chance that Q{w) is equal to any particular 
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standard Young tableau of shape A is sa • • • , jf) , where there are b^' variables. Letting 
/a (a) denote the number of standard Young tableaux of shape A with a descents, it follows 
that 

\\\=n ^ ^ 

By equation 7.96 of [25\ . 



fxiar^' = (1 - tr^' Yl • • - 1)*' 

a>0 k>l 

where in the kth summand, sa(1) • • • , 1) denotes the Schur function with k variables spe- 
cialized to 1. Thus 



^P(dH = a) •t'^+i 

a>0 

a>0\X\=n ^ ^ 

= {i-tr'Y^'Y^^ .a(i,...,i) 

fe>l |A|=n ^ ^ 



A;>1 n>0|A|=n 

where [u"']g{u) denotes the coefficient of is a power series g{u). Applying the Cauchy 
identity for Schur functions |25l p. 322], this becomes 

k>l k>l V / 

as desired. □ 
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